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Spin glass order in spin turbulence 
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We study numerically the spin turbulence in spin-1 ferromagnetic spinor Bose-Einstein conden- 
sates (BECs). Spin turbulence (ST) is characterized by a —7/3 power law in the spectrum of the 
spin-dependent interaction energy The direction of the spin density vector is spatially disordered 
but temporally frozen in ST, showing some analogy with the spin glass state. Thus, we introduce the 
order parameter of spin glass into ST in spinor BECs. When ST develops through some instability, 
the order parameter grows with a —7/3 power law, thus succeeding in describing ST well. 
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Most flow in nature is actually turbulent [1,]. As the 
velocity increases, the flow generally changes from lam- 
inar to turbulent. Turbulence has been a great mystery 
for a long time mainly because it is a complicated dy- 
namical phenomenon with strong nonlinearity far from 
an equilibrium state. When we consider such compli- 
cated phenomena, it is important to focus on statistical 
laws. The most important statistical law in turbulence 
is the Kolmogorov —5/3 law. The Kolmogorov energy 
spectrum refers to E{k) = Ce^/^/c~^/^, where k is the 
wave number, e is the energy flux, and C is the Kol- 
mogorov constant of the order unity. The Kolmogorov 
spectrum is confirmed in fully developed turbulence in 
classical fluids, being the smoking gun of turbulence. 

Apart from these studies in classical fluid dynamics, 
there have been studies on quantum fluids such as su- 
perfluid helium and atomic Bose-Einstein condensates 
(BECs). These systems are subject to quantum restric- 
tions; the appearance of order parameters means that 
rotational motion can exist only in the presence of quan- 
tized vortices. A quantized vortex is a stable topological 
defect with quantized circulation. Such vortices give rise 
to quantum turbulence (QT) [2|, |3[. Because quantized 
vortices are well defined as elements composing a turbu- 
lent flow, QT is expected to be an easier system to study 
than classical turbulence and gives a simpler model of 
turbulence. Quantum turbulence has long been studied 
in superfluid helium, whereas atomic BEC has only re- 
cently become an important research area. 

An important interest on QT is how to characterize 
the turbulent state. Numerical simulation of the Gross- 
Pitaevskii (GP) model shows that the spectrum of im- 
compressible kinetic energy obeys the Kolmogorov —5/3 
law in a uniform system [j, |5| and in a trapped system 
[6[ . Confirmation of such an energy spectrum power law 
would certainly provide strong proof of turbulence. How- 
ever, there must be another better way to characterize 
QT. 

Spinor BECs have spin degrees of freedom and exhibit 
phenomena characteristic of spin [7] . The hydrodynamics 



of spinor BECs has recently been studied by several au- 
thors [sMl 1| . In previous studies, we found that hydrody- 
namic instability in spin-1 spinor BECs exhibits unique 
behavior owing to their spin degrees of freedom and form 
spin turbulence (ST) in which the spin density vector 
has various directions. Although the spin-dependent in- 
teraction can be ferromagnetic or antiferromagnetic, we 
confine ourselves to the case of ferromagnetic interaction 
in this work. The first study addressed the counterflow 
between the m = 1 and m = —1 components in a uni- 
form system, where m is the magnetic quantum number. 
Through the instability the counterflow leads to ST, in 
which the spectrum of the spin-dependent interaction en- 
ergy obeys a —7/3 power law [12]; the —7/3 power law is a 
statistical law characteristic of ST and can be understood 
by the scaling analysis of the time-development equation 
of the spin density vector. Such a spin-disordered state 
was experimentally created in a trapped system through 
the instability of the initial helical structure of spins [13| , 
which allowed us to study ST numerically in a similar 
situation and confirm the —7/3 power law again [14]. An 
oscillating magnetic field applied to a uniform ferromag- 
netic system can also create ST with a —7/3 power law 

These three works 12|, [ij, ll5| reveal the important 
characteristics of ST. First, the —7/3 power law is ro- 
bust independently of whether the system is uniform or 
trapped. How the system is excited or how the ST is 
created has little effect on the power law. Second, ob- 
servation of the spin motion in all three cases indicates 
that the spin density vectors become spatially random 
but temporally frozen as the ST develops and the —7/3 
power law appears clearly, which reminds us of spin glass. 
Spin glasses are magnetic systems in which the interac- 
tions between the magnetic moments are in conflict with 
each other [l6|. Thus, no conventional long-range order 
can be established. Nevertheless, these systems exhibit 
a freezing transition to a state with a kind of order in 
which the spins are aligned in random directions. 

In this paper, we introduce the order parameter of spin 




FIG. 1: (Color online) Distribution of spin density vectors in 
ST at tlr — 4000 obtained by the counterflow instability with 
yRJCsound = 0.78, where Vr, Csound-, and r are the relative 
velocity of counterflow, the sound velocity, and the charac- 
teristic time in a uniform system [121 . [2l| | . The system size is 
128^ X 128^, where ^ is the coherence length. This numerical 
calculation starts from the initial state of the counterflow be- 
tween the vfi—\ and —1 components, which induces the spin 
modulation, finally leading to ST. Details of the dynamics are 
described in \\%. 



glass [16, 17] to characterize ST in spinor BECs. This 
introduction is so successful that the order parameter is 
found to grow with the setup of the —7/3 power law. This 
success paves the way for two innovative approaches. One 
is to propose a useful order parameter in the turbulence 
of quantum fluids; such a proposal is impossible in ordi- 
nary QT [2| . The other is to connect the study of spinor 
BECs with that of magnetism including spin glass. In 
this paper, first we describe how the ST in spinor BECs 
behaves. As the directions of spins become disordered, 
the spectrum of the spin-dependent interaction energy 
starts to exhibit the —7/3 power law and the motion of 
the spins becomes temporally frozen. Second, the intro- 
duction of the spin glass order parameter is shown to 
successfully characterize the ST. 

We consider a spin-1 spinor BEC at zero temperature, 
which is well described by macroscopic wave functions 
'0m {j^ = 1, , — 1). The wave functions '^rn obey the GP 
equation [11, [l9| 
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n=-l 

1 
+ ComlJm + Ci ^ S • Smn'4^n- (1) 

n=-l 

Here, V and B are the trapping potential and mag- 
netic field. The parameters M, ^, /i^, and q are the 
mass of a particle, the Lande g factor, the Bohr mag- 
neton, and a coefficient of the quadratic Zeeman effect, 
respectively. The total density n and the spin density 
vector Si {i = x,y,z ) are given by n = Z)m=-i IV^^P 
and Si = YlL,n=-i'^mi^i)mn'^n with the spin-1 matri- 



ces {Si)mn- The parameters cq and ci are the coefficients 
of the spin-independent and spin-dependent interactions. 
In our research of ST, we focus on the spin-dependent in- 
teraction energy Egpi^ = ^ J s'^dr^ whose coefficient ci 
determines whether the system is ferromagnetic (ci < 0) 
or antiferromagnetic (ci > 0). We are interested in the 
ferromagnetic case where ST clearly exhibits the —7/3 
power law. 

In this paper, we study the two-dimensional ST ob- 
tained by (i) the counterflow instability [12], (ii) the in- 
stability of spin helical structure [14], and (iii) an oscillat- 
ing magnetic field [l5| . All parameters of the numerical 
calculation of (i) and (ii) in this paper are the same as in 
the previous studies. In the case of (iii), we use parame- 
ters different from those of the previous study [20| 

Figure 1 shows a typical case of ST developing from the 
instability of the counterflow between the m = 1 and m = 
— 1 components [12]. As ST develops, the spectrum of the 
spin-dependent interaction energy begins to exhibit the 
—7/3 power law. Figure 2 shows that the —7/3 power law 
appears in all three cases we have studied: counterflow 
instability in a uniform system [Fig. 2(a)] [12], instability 
of the initial helical structure of the spin density vector 
in a trapped system [Fig. 2(b)] [14], and application of 
an oscillating magnetic field in a uniform ferromagnetic 
system [Fig. 2(c)] [15]. Thus, the —7/3 power law is 
universal in ST independently of whether the system is 
uniform or trapped or how the ST is created. All cases 
show that the spin density vector is spatially random but 
temporally frozen; the case of the counterflow instability 
is shown in the movie [21]. We know such a phenomenon, 
namely spin glass, which invites us to introduce the order 
parameter of spin glass. 

The order parameter of spin glass is ordinarily intro- 
duced as follows: An equilibrium system is supposed 
to consist of TV lattice sites, each site i having spin Si 
{i = 1, A/"). The time average {Si{t)) and the space aver- 
age [Si{t)] are defined as 



{sm 



lim — 



^ Si{t)dt, [S,(t)] = l^Si(t). 



(2) 
Then it is possible to introduce two order parameters, 
namely, the mag netization M = [{Si{t))] and q = 
[{Si{t))'^] [16|, ll7|. If the system is paramagnetic, both 
M and q vanish. A ferromagnetic order gives nonzero M 
and q. When M = but q is nonzero, the state is called 
spin glass [17], and the directions of spin are spatially 
random but temporally frozen. 

To apply the order parameter to our present case, 
we have to consider two things. First, we address the 
spin density vector s(r,t) instead of Si{t). Our spinor 
BECs are usually trapped so that the amplitude of the 
spin density vector is not uniform. To extract this ef- 
fect and focus on the direction of the spin density vec- 
tor, we define the order parameter with the unit vector 
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FIG. 2: (Color online) Spectrum of the spin-dependent interaction energy in three cases: (a) counterflow instability in a uniform 
system, (b) instability of the initial helical structure of the spin density vector in a trapped system, and (c) application of an 
oscillating magnetic field in a uniform ferromagnetic system. Each spectrum is calculated at (a) t/r = 4000, (b) tu = 200, 
and (c) t/r = 4000, when ST is fully developed. The parameters uj, an, Rtf, no, ^, L, and Am are, respectively, the trapping 
frequency, the harmonic oscillator length, the Thomas-Fermi radius, the total density in the uniform system, the coherence 
length, the system size, and the wavelength where the energy is injected by the oscillating magnetic field. The details of these 
parameters except for Am are described in our previous paper [12. \14, \1^ . The expression for Am is obtained by the resonant 
condition with the oscillating magnetic field [22]. The blue broken shows the —7/3 power law. The spectrum in all cases 
exhibits the —7/3 power law in the scaling region [(a) 27v/L < k < 2n /^s , (b) 2'k/Rtf < k < 27r/.fs, (c) 27v/Xm < k < 27r/^s]. 



s(r,t) = s(r,t)/|s(r,t)|. The definition of the space av- 
erage of Eq. (2) is replaced by 



[s{r.t)] = jl s{r,t)dr 



(3) 



with system area A. Second, the system starts from some 
initial state to develop ST, which should be characterized 
by the time-dependent order parameter. Thus, we intro- 
duce the time average of s(r, t) during the period [t, t+T] 
as 



{s{r 



1 f'^^ 



s{r^ti)dti 



and define the time-dependent order parameter 
q{t) = [{s{r,t)fr]- 



(4) 



(5) 



When we calculate q{t) for our system, we should be 
careful with how we take T. Generally, a longer T is 
desirable, but in reality, we should estimate q{t) with 
some finite T. The criterion for the appropriate value of 
T would be that T should be longer than some charac- 
teristic time of the system. In this system, the velocity 
characteristic of spin is given by Cg = y^\ci\n/M [23]. 
Then, a typical characteristic time is the system size L 
divided by the velocity Cg . We will show in the following 
that this time L/cg is found to be so long that the system 
in the initial state becomes sufficiently disturbed by the 
time. Thus, we take T comparable to L/cg. In ST, the 
order parameter q{t) decreases with T, but the depen- 
dence is weak. If the spin density vector is completely 
frozen in ST, q{t) should be unity. 

We show how q{t) grows toward ST with a —7/3 power 
law of the spin-dependent interaction energy. Figure 3 



shows the time dependence of the power exponent of the 
spectrum of the spin-dependent interaction energy ob- 
tained by the least-squares method and q{t) for (a) coun- 
terflow instability in a uniform system, (b) instability of 
the helical structure in a trapped system, and (c) appli- 
cation of the oscillating magnetic field in a uniform fer- 
romagnetic system. Figures 3(a) and 3(b) show that q{t) 
increases obviously as the exponent n approaches —7/3. 
Furthermore, the magnetization \m{t)\ = |[(s(r, t))T]| is 
much smaller than \/q{t) because [s(r, t)] is conserved to 
keep vanishing under Eq. (1) without the magnetic field. 
Thus, the two order parameters q{t) and \m{t)\ are found 
to be effective for describing the ST in the two cases. 

In contrast, q{t) of Fig. 3(c) behaves differently from 
that of Figs. 3(a) and 3(b). In the case of Fig. 3(c), 
[s(r,t)] is not conserved; in Fig. 3(d) we see that \m{t)\ 
decreases from the initial ferromagnetic state toward 
ST. In ST, ^/q(f) of Fig. 3(c) is smaller than that of 
Figs. 3(a) and 3(b), because the oscillating magnetic 
field causes fluctuation in the spin density vector and re- 
duces \/q{t). However, y/q{t) ~ 0.35 is much larger than 
the magnetization \m{t)\ ~ 0.005 at t/r = 4000, which 
means the growth of the spin glass order in ST. There- 
fore, the two order parameters q{t) and \m{t)\ are still 
found to be effective for describing the ST in this case. 

Only case (b) is affected by the inhomogeneity. Be- 
cause the system is trapped, the amplitude of the spin 
density vector is reduced going from the center to the 
boundary. The low condensate density near the bound- 
ary induces a large fluctuation of the condensate phase, 
thus causing s(r,t) to fluctuate temporally. As a result, 
the space average of Eq. (|3|) depends on the area A of 
the integral. Figure 3(b) shows two cases for which the 




1 

0.8 
0.6 



0.4 \y^ 



- 0.2 



5000 



(b) 




1 




4 


0.8 




2 



0.6 


n 


-2 
-4 


0.4 




-6 


0.2 




-8 
■10 



-12 



150 200 



i ' ' 0.1 -♦ 

-I 500 1000 


-: n = -7/3 



2500 



1 

0.8 

0.6 

0.4 

0.2 





V^ 



5000 



FIG. 3: (Color online) Time development of the exponent n, the spin glass order parameter y^q{t), and \m{t)\. n and \/q{t) 
are shown for (a) the counterflow instability in a uniform system, (b) the instability of the initial helical structure of the spin 
density vector in a trapped system, and (c) the application of the oscillating magnetic field in a uniform ferromagnetic system. 
The time development of |m(t)| in the case of (c) is shown in (d). In (b), the squares and triangles show that the calculation 
of -yjq is performed in the regions r < Rtf /'^ and r < Rtf^ respectively. The exponents n in (a) ^ (c) are evaluated by the 
least-squares method in the scaling region defined in the caption of Fig. 2. 



radius of A is Rtf/'^ and Rtf with the Thomas- Fermi 
radius being Rtf- The larger value of q{t) for Rtf l'^ 
means that the central spins are more likely to be frozen. 

Finally, we note that the time T = Lj Cs is so long that 
the system in the initial state becomes disturbed, which 
is found by the time development of n. The time T is 
(a) 572r, (b) 24/cj, and (c) 572r, respectively. We find 
that n rapidly changes in the period < t < T, which 
means that the distribution of the spin density vector 
in the wave number space changes rapidly too. Thus, in 
this period, the spin density vector temporally points in 
various directions. This is also confirmed in the movie 
of the spin density vector [21]. Therefore, the time L/cg 
is appropriate for the calculation of the spin glass order 
parameter q(i). 

In summary, we have studied numerically ST in spin-1 
spinor BECs. The spectrum of the spin-dependent in- 
teraction energy is found to exhibit a —7/3 power law 
independently of the details of the system or how the 
ST is created. The direction of the spin density vec- 
tor is spatially disordered but temporally frozen in ST, 
which reminds us of the spin glass state. Following spin 
glass theory, we introduced the order parameter of the 
spin glass into ST in spinor BECs. The order parameter 
q(i) of spin glass grows with a —7/3 power law, well de- 
scribing ST. Therefore, ST can be characterized by three 
quantities: (i) the spin glass order parameter yjq^ (ii) the 
magnetization m, and (iii) the power exponent n = —7/3 
of the spectrum of the spin-dependent interaction energy. 
These behaviors should be accessible experimentally. 

We thank Shin-ichi Sasa for useful discussions. 
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